Critical current under an optimal time-dependent polarization direction for Stoner 
particles in spin-transfer torque induced fast magnetization reversal 
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Fast magnetization reversal of uniaxial Stoner particles by spin-transfer torque due to the spin- 
polarized electric current is investigated. It is found that a current with a properly designed time- 
dependent polarization direction can dramatically reduce the critical current density required to 
reverse a magnetization. Under the condition that the magnitude and the polarization degree of 
the current do not vary with time, the shape of the optimal time-dependent polarization direction 
is obtained such that the magnetization reversal is the fastest. 
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The advent of miniaturization and fabrication of mag- 
netic particles at nano-meter scale (called Stonerparti- 
cles) makes the Stoner- Wohlfarth (SW) problem|3 very 
relevant to nano-sciences and nano-technologies. For a 
Stoner particle, strong exchange interactions keep the 
magnetic moments of atoms rigid such that the con- 
stituent spins rotate in unison. The magnetization 
dynamics can be manipulated by a laser light 3] , or 
a magnetic field lj , M , or a spin-polarized electric 
current [111 II [iflEMlEllillii- Among them, mag- 
netization manipulation 13 of Stoner particles by electric 
current of spin-polarized electrons is of great current in- 
terests in spintronics because of its locality and low power 
consumption. The idea of spin-transfer torque (STT) 
generated by a spin-polarized electric current was in- 
dependently suggested[3 by Slonczewski and Berger in 
1996, and was verified by several experiments [lOj . Im- 
portant issues in its applications are to lower the critical 
current required to reverse a magnetization^^ and to 
design a current pulse such that the magnetization can 
be switched from one state to another fast. Many rever- 
sal schemes IH E E |Il|l3 have been proposed and 
examined. Ideas include thermal- assistance ^3 E| and 
sample designs [isl lisj . However, the question of how 
much the critical current can further be lowered is still 
unknown. In this report, the optimal time-dependence of 
the current polarization direction for the fast magnetiza- 
tion reversal is derived when the magnitude and the po- 
larization degree of the current is fixed. It is shown that 
the critical current density can be dramatically reduced 
if the optimal time-dependent polarization is employed. 

The prototype of STT systems is a magnetic multi- 
layered structure of nano-meter scale as illustrated in 
Fig. 1. It consists of two ferromagnets which are sand- 
wiched among three nonmagnetic metallic layers. Elec- 
trons travel along the +z direction (from the left to the 
right) in the sample (the direction of the current is op- 
posite as shown in Fig.l). The first ferromagnet Fl is 
usually very thick so that the current does not affect the 
magnetization Mi of Fl. Electrons are polarized along 
Ml after they pass through Fl, and remain their polar- 
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FIG. 1: Schematic illustration of the STT structure. Note 
that the direction of the electrical current is opposite to that 
of electron flow. 



izations before entering the second ferromagnet F when 
the thickness of the spacer layer between Fl and F is 
much smaller than the spin diffusion length. The polar- 
ized electrons transfer their spin angular momenta to F, 
resulting in so-called STtQ. This STT can affect the 
dynamics of magnetization M of F when its thickness is 
thin enough. Theoretical studies^, S show that the 
STT r is proportional to the current with the following 
form 



^ = ISTT = g{P,m- s)mx {mx s) (1) 

at fioe 



where m, s are the unit vectors of M and the polariza- 
tion direction (along Mi) of the current, respectively. 
In the expression, V, h, fiQ, and e denote the volume 
of F, the Planck constant, the vacuum magnetic per- 
meability, and the electron charge, respectively. 7 = 
2.21 X 10^ {rad/ s)/ {A/ m) is the gyromagnetic ratio. The 
exact microscopic formulation of the STT is still a sub- 
ject of study and debateH^]. Different theories differ 
themselves in different expressions of function g that de- 
pends on the degree of the polarization P of the current 
and relative angle between m and s. All experimental 
investigations 10] so far are consistent with the result of 
Slonczewski ,7] which will be used throughout this study. 
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9{P. • ^) = (1 ^ p)3(3 + rn-s)- 16P3/2 " 

The magnetization dynamics of a Stoner particle (F 
in Fig. 1) under an effective magnetic field Ht and a 
polarized current is governed by the so-called modified 
Landau-Lifshitz- Gilbert (LLG) equationjsj with an addi- 
tional term (the third term on the right of the equation 
below) due to the STT of Eq. ((TJ 

-— = -7Af X Ht+amx — — + jaiM x {M x s), (3) 

where a/ = hig/ {fiQcAPV) is a dimensionless parame- 
ter, and a is a phenomenological dimensionless damping 
constant whose typical value ranges from 0.01 to 0.22 
for Co films Because the magnitude of M does not 
change with time for a Stoner particle, M can be de- 
scribed by the polar angle 9 and the azimuthal angle 
in the spherical coordinate, and Eq. ^ can be rewritten 
in a dimensionless form 

(1 -I- a^)^^ ~ —rh X hi ~ m X (m X /12), (4) 
at 

where 

hi ^ ht + aajs, 
/i2 — OLht — ais. 

t in Eq. |0J is measured in unit of {jA4)^^. The mag- 
netization and the magnetic field are in the units of M. 
the total field ht = h + hi includes both the applied mag- 
netic field h and the internal field hi due to the magnetic 
anisotropic energy density ■u;(m), hi — ~V^w{rh)/ ^q. 
Different particles are characterized by different magnetic 
anisotropy. We consider only the uniaxial particles of 
magnetic anisotropy w{rh) — — fcm^/2 with its easy axis 
along the z-direction so that and hi ~ ^^E^EE^Wy^ = 
f{cos0)z. Let Cr, ee, be the three spherical unit vec- 
tors of m, so hi = — fcsinflcosflee -|- fccos^ 9er. In terms 
of 9 and (p, Eq. Q can be written as 

(1 + a'^)9 ht^4, + ahtfi + ai{as^ - sg), 

(1 + a^) sa\9(t) — ahf — ht^g — ai{asg + s^). 

Here ht.e, ht.4, and sg, s^, Sr are the eg, icf,, and ir com- 
ponents of ht and s, respectively. 

Because hi and /12 are in general non-coUinear, the dy- 
namics with the additional STT term in Eq. 10} is very 
different from that without this term which describes a 
Stoner particle only in a magnetic field. The particle en- 
ergy can only decrease in a static magnetic field since 
the field cannot be an energy source0. However, a po- 
larized electric current can pump energy into the particle 



through the STT. Thus, the STT term from even a dc 
current can be an energy source, and the dynamical be- 
havior of a Stoner particle under an STT is much richer 
than its counterpart in a static magnetic field. According 
to Eq. the magnetization undergoes a precessional 
motion around field hi and a damping motion toward a 
different field /i2- 

The magnetization switching problem for a uniaxial 
Stoner particle is as follows: In the absence of both polar- 
ized current and h{— 0), M in F has two stable directions, 
9 — and 9 — tt, along its easy axis (z-axis). The goal is 
to reverse the initial state (say 9 = 0) to the target state, 
9 — TT, fast by a small polarized electric current. All stud- 
ies so far assumed that the current polarization direction 
does not vary with time. However, previous studies |^ on 
magnetic-field induced magnetization reversal show that 
a switching field can be dramatically reduced if the direc- 
tion of the field varies properly with time during a rever- 
sal. Thus, it is natural to investigate whether one can use 
a current with a proper time-dependent polarization di- 
rection to lower the critical reversal current density. The 
precise issue studied here is: For a given system in which 
damping constant a, anisotropy f{cos9) and the exter- 
nal magnetic field are fixed, the polarized electric current 
can vary its polarizations direction during a reversal pro- 
cess under the constrain of the constant magnitude / and 
constant polarization degree P of the current. The issue 
is to minimize the critical current density, and to find 
the best shape of polarization direction for the shortest 
reversal time. 

To simplify the calculations, we consider the case of 
zero external magnetic field, h = 0. However, the basic 
idea should directly be applicable for non-zero static ex- 
ternal field because h can be added to hi. Then Eq. © 
becomes, 

(1 + a'^)9 = ai{as^ — sg) — a/(cos 9) sin 9, 

■ . (6) 

{1 + a ) sin 9(f) = —aj{asg + s^) + /(cos 9) sin 9. 

s — s{t) is in general a function of the time. Differ- 
ent function will lead to different angular velocities for 
9 and 0. Thus, the magnetization reversal time from 
the initial state {9 — 0) to the target state {9 — tt), 
defined as T = d9/9, depends on s{t). In order to 
find the optimal s{t) that minimizes T, one only needs 
ai{as^ — Sg) or g{P, Sr){as^ — sg) to be maximum such 
that 9, according to Eq. will be the largest at any 9. 
This observation is important and it can be applied to 
other function forms of g. Because + + s"^ = 1, the 
maximum of g{P, Sr)(Q!S0 — sg) can be obtained from the 
standard Lagrange multiplier method in which one intro- 
duces L = g{P, Sr){as^ — sg) — X{sl + Sg + s'^). By setting 
the partial derivatives of L with respect to Si (i=r, 9, (p) 
to zeros, the maximum of g{P, Sr){as,f, — sg) is 

[g{P, Sr){as^ - se)]max = VTT^G(F), (7) 
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where G(P) = g{P, s*)^^! - sf and 



16P3/2 _ 3(1 + p)3 
1 



a r- 



*2 
r ' 



(8) 



Eq. ijHJl gives the optimal polarization direction which will 
lead to the shortest switching time under a fixed current 
magnitude. Although the optimal s* appears to depend 
only on damping constant a and P, but not on / (cos 6), it 
is in fact time-dependent since Cr, eg, vary with the 
time. Furthermore, magnetic anisotropy /(cos 6*) shall 
influence the evolution of rh which in turn influences the 
time-dependence of s*. Thus, if they were to change 
/(cos 6*) and nothing else, the time-dependent s* would 
be different. 

Under the optimal design of Eq. ©, 9{t) and (f>{t) will 
satisfy, respectively. 



hi 



G{P) 



af{cose)sm0/{l + a^), (9) 



and 



= /(cos0)/(l + a2). 



(10) 



For w{'m) = — it is straightforward to integrate 

Eq. and obtain the reversal time T from 6 = to 

6* = TT, 



T - 



{o? + l)7r 



(11) 



In the weak damping limit (a — > 0) or large current limit 

(/ OO), T (X Tl/ 1. 

For a uniaxial model, the critical switching current or 
current density Jc can be obtained by setting 6 — Q in 
Eq. (jni . This is because 9 cannot be negative if the mag- 
netization of a uniaxial particle moves from 9 — Q to 
9 — TT. Therefore, the first term in Eq. ^ must ex- 
ceed the second term due to magnetic anisotropy for all 
9's {g [0, tt]) in a reversal. This simplicity for a uniax- 
ial model comes from ^-independence of Eq. I^. Since 
Eq. is the largest velocity for an arbitrary 9 under 
the best choice of the polarization direction of a current, 
the critical current should be the one when the smallest 
9 (for all 9) is zero. The critical reversal current density 
in our case is 



Jr = 



hG{P) VTT^ 



Q. 



(12) 



Here the current density is defined as J = //A with A 
being the cross-section area and d being the thickness 
of F. Q = maj:{/(cos6')sin6l} for 9 e [0,ti]. Q = k/2 



for w{rri) = —krnl/2. One notices, from the derivation 
of Eq. I|12|) . that the assumptions of time- independence 
of / and P are not essential for Jc as long as STT is 
proportional to / and g in Eq. (1). 

To see how much the critical switching current density 
is reduced, let us compare the critical current density 
of Eq. H12|l with those in the previous schemes where 
the current polarization direction is fixed. When the po- 
larization direction s is parallel to the easy axis of the 
uniaxial magnet F (the critical current density in the 
perpendicular case is larger), the critical current density 
for the same anisotropy as that of Eq. (|12|l is0,0,0| 



Jr. = 



hgiP, 1) ' 



(13) 



Fig. 2 is the plot of Jc versus damping constant a for 
P = 0.4. The dashed line is that for Eq. and the 

solid line is the result of our new strategy which saturates 
to a constant at large a limit. At a = 0.1, Jc in the new 
strategy is about one fourth of that given by Eq. H13|l . 
The difference between Eqs. (|12f) and H13I) depends on 
the degree of polarization, P. Fig. 3 is Jc vs. P at a = 0.1. 
It should be pointed out that zero Jc in Eq. H12|l at P = 1 
is an artificial result originated from the divergence of 
g{l,x) at X = —1 in the Slonczewski's theorylTl. This 
divergence is removed in other formulations of g 9\. 




FIG. 2: Jc vs. a for P = 0.4 and uniaxial model of w{m) = 

Although our results are obtained for uniaxial Stoner 
particles with a specific function of g, the basic ideas 
and approaches can in principle be generalized to the 
non-uniaxial cases with other ^-functions. Also, almost 
all experiments so far used a set-up like that illustrated 
in Fig.l. However, the formulation of an STT induced 
magnetization reversal does not really rely on the set- 
up. In general, the same formulation could be applied 
to a Stoner particle through which a spin current pass. 
Therefore, advance in generating the pure spin current 
(with charge current being zero), one of current topics 
in so-called spin Hall effect, could lead to a new set of 
experimental set-up beyond the prototype illustrated in 
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FIG. 3: Jc vs. P at a = 0.1. The rest of system parameters 
are the same as that of Fig. 2. 



Fig. 1 in the STT experiments and the STT apphca- 
tions. In the prototype of Fig. 1, it happens that the 
current of spin-polarized electrons is generated by one of 
the magnets. 

The realization of the results reported here depends 
largely on our ability in generating a designed spin- 
polarized electric current. In some sense, the present 
work converts the issue of critical current to the issue of 
generating an arbitrary polarized electric current. Any 
breakthrough in the front of spin current generation shall 
lead to the great leap forward in magnetization minipu- 
lation. One possible way to generate a desired polarized 
electric current is through controlling magnetization Mi 
of Fl in Fig. 1 by other means. As it was explained 
early, the polarization direction of the electric current is 
parallel to Mi in an experimental set-up illustrated in 
Fig. 1. Thus, time-dependence of the polarization of the 
current is the same as that of Mi. However, since the 
response time of electrons to an electric current is usu- 
ally much faster than that of the magnetization, it may 
be an experimental challenge to create a required time- 
dependent spin polarization by the method, especially 
when the change of polarization is very fast. 

In conclusion, we have showed that a proper time- 
dependent polarized electric current can dramatically re- 
duce the critical current density needed to reverse a mag- 
netization. An optimal time-dependent current polariza- 
tion is obtained such that the magnetization reversal time 
is the shortest for uniaxial Stoner particles. 
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